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ABSTRACT 

We generalize in several ways the results existing in the literature: a) we make 
use of an exact general relativistic solution for a spherical, nearly empty cavity 
in the matter dominated era to evaluate the null geodesies and the Sachs- 
Wolfe effect; b) we evaluate the magnitude of the adiabatic fluctuations of the 
photon-baryon plasma; c) we study the influence of the shell profile; and d) we 
take into account the finite thickness of the last scattering surface (LSS) and 
the influence of its position with respect to the void center. We find empirically 
an analytic approximation to the Sachs- Wolfe effect for all crossing geometries 
and we derive an upper limit of w 25 h~ l Mpc for the comoving radii of voids 
sitting on the LSS in order to achieve compatibility with COBE's data. As a 
nearly empty void has an overcomoving expansion of a factor of rs 4 between 
decoupling and the present, the maximum allowed size at present is ks 100 
h^ 1 Mpc. On the other hand, the smallness of the comoving size relative to 
the sound horizon reduces strongly the adiabatic effect by Silk damping and 
makes it negligible. Most of the signature of primordial voids comes therefore 
from metric effects and consists of sub-degree spots blue or red depending on 
whether the center lies beyond or within the LSS. In conclusion we refine 
and confirm earlier constraints on a power law void spectrum originated in 
an inflationary phase transition and capable of generating the observed large 
scale structure. 



INTRODUCTION 



Inflation provides two well known mechanisms for relat- 
ing primordial physics to perturbations on large scales 
(> 10 /i _1 Mpc) in the observed universe. The first one 
arises from the small quantum fluctuations of the in- 
flaton field: the fluctuations are stretched out to scales 
greater than the horizon (H^ 1 ) by the spatial expan- 
sion, and their dynamics freezes until the reenter in 
post-inflationary epochs. The second mechanism oc- 
curs in the case of a first order phase transition in the 
early universe; the nucleated bubbles imprint the den- 
sity field. In fact, in this work we assume that they are 
the primordial cause of the large scale voids observed 
at the present. A class of inflationary models provides 
bubbles production; the underlying minimal physics is 
given by quantum field theory through quantum transi- 
tion of the inflaton field from a false vacuum state to a 



true vacuum one. More realistic models of first order in- 
flation (FOI), after the original version by Guth (1981), 
were found by considering non-minimal gravity theories 
and two-field phenomenology (see e.g. La & Steinhardt 
1989, Kolb 1991, Adams & Freese 1991, Occhionero & 
Amendola 1994), and made FOI the richest extension 
of the slow-rolling model worked out by Linde (1983). 

The hypothesis of a bubbly distribution of mat- 
ter in the universe, theoretically interesting in itself, is 
strongly suggested by recent observations. Large-scale 
voids (> 10 /i -1 Mpc ) have been discovered by recent 
redshift surveys (see e.g. Kauffmann & Fairall 1991, Vo- 
geley, Geller & Huchra 1991, Lindner et al. 1995, and 
references therein), but the definition of a statistics for 
their dimensions is still premature, due to the pres- 
ence of voids comparable with the sample's size. Deep 
(1000 h~ Mpc ) pencil-beam surveys toward the galac- 
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tic poles revealed a fascinating regularity in the galaxy 
distribution with a characteristic scale of 128 /i _1 Mpc 
(Broadhurst et al. 1990); El- Ad et al. (1996) constructed 
an algorithm to find voids in redshift surveys; da Costa 
et al. (1996) found that the voids are really empty of 
matter. 

The estimate of the anisotropies induced by a void 
in the cosmic microwave background (CMB) yields 
the strongest constraint to the underlying inflationary 
physics by comparison with COBE observations (see e.g. 
Bennet et al. 1994 and references therein). Detections 
of degree scale CMB anisotropies (see e.g. Scott et al. 
1996) provide in our case weaker constraints because 
they involve a limited part of the sky. To be quanti- 
tative, upgrading the experiment from 10 to 1 degree 
would imply a gain of a factor 100 in resolution; but 
in our case we would be looking for one single object, 
the largest, and we have to take into account the por- 
tion of sky in which we search (precisely rescaling by 
the corrisponding factor cLNb in Eq.(37) below). Conse- 
quently, such observations would be significative in our 
case if the covered solid angle were at least 400 square 
degrees, and to our knowledge high resolution observa- 
tions of this width do not exist. 

The first attempt to relate the evidence of voids in 
the large scale matter distribution with FOI was made 
by La in 1991. Overcoming the difficulties of previous 
works (Turner et al. 1992, Liddle & Wands 1991, 1992), 
it has been found that it is possible to construct mod- 
els that fit the galaxy correlation function and pass the 
CMB bounds with a bubbly matter distribution (Amen- 
dola & Occhionero 1993, Amendola & Borgani 1994). 
Occhionero & Amendola (1994) found a toy model of 
FOI where the expected distribution in radius of the 
primordial voids, defined as the number of voids with 
radius greater than R today, is approximated by a power 
law, 

N B {R)=(^f) P . (1) 

This distribution is compatible with the measured 
galaxy power spectrum and with the CMB anisotropies 
for 

R M ~ 30 /i _1 Mpc , 6 < p < 10 , (2) 

normalizing Nb to a cube with 500 h~ Mpc by side. 
Aim of this paper is to review and to set the constraint 



on Rm and p more rigorously than in the original paper, 
by an accurate evaluation of the physics involved. 

A void is characterized by its central to asymptot- 
ical density ratio 

^ Pcenter 

Poo 

that decreases with time. The dynamics of voids has 
been investigated by Icke (1984) and Bertschinger (1983, 
1985) in the strongly non-linear regime, A 1, and for 
t — » co. The important result of their analysis is that the 
physical radius r of a compensated void of total energy 
E increases as £ 4 / 5 and that all the significative quan- 
tities of a void (pressure, density profile, velocity field, 
mass) are only functions of the dimensionless coordi- 
nate A = r(Et 2 1 Poo) -1 '' 5 , thus remaining unchanged at 
late times (self-similar expansion). These analyses have 
been confirmed by iV-body simulations (White & Os- 
triker 1990, Dubinski et al. 1993). In the following we 
mean that R in Eq. (^) is the physical radius today; how- 
ever, the corrisponding comoving quantity was smaller 
in the past by a factor of (t/t )^~i = (1 + z) -0 ' 2 ; in 
particular, a void on the last scattering surface (LSS), 
at 1 + z ~ 1000, was physically smaller by a factor of 
4 x 1000. 

If we connect present voids with primordial bub- 
bles, the energy profile at nucleation time can in princi- 
ple be calculated exactly by solving the inhomogeneous 
Klein-Gordon equation for the inflaton field with the 
appropriate tunneling-like initial conditions. The gen- 
eral structure is a spherically symmetric underdensity, 
with a thin wall, and a size of the order of H -1 (Cole- 
man 1977). The internal matter content depends on the 
detailed shape of the transition potential. However, its 
successive evolution depends crucially on the proper- 
ties of the dominating cosmic fluid. From Vadas (1993) 
we know that super-horizon voids, during the radiation 
dominated era (RDE), undergo hydrodynamical inflow 
of the primordial plasma into the empty central region. 
For an estimate of this effect we call r/ the mean veloc- 
ity dispersion of the primordial plasma in units c = 1, 
where the cases rj — 0,1 mean respectively predomi- 
nance and absence of CDM; simply, we sketch the fill- 
ing process as inflow of cosmic fluid in a spherical empty 
void of comoving size R/4 (the 4 factor accounting for 
the overcomoving expansion) from the end of inflation 
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(T e ~ 10 15 GeV) across equivalence (Teq — lOeV) un- 
til decoupling (T D ~ 2 ■ 10 _1 eV); one finds that the 
quantity A defined in (H) is given by 



A(fo) 



exp 



-v 



2H~ 



R/4 V t d 



T E q ) 



•(4) 



For R ~ 10 /i 1 Mpc our estimate gives A ~ 1 for 77 = 
O(ijjo), a very typical value for CDM (La 1991). Thus, 
voids smaller than J? ~ 10 /i _1 Mpc are erased by matter 
inflow; this is therefore a lower cut-off in the distribution 
(0). We can see that larger voids (R ~ 10^50 /i _1 Mpc) 
have a typical filling fraction .5 > A > .1 in their central 
underdense region. Since no significative dependence of 
our results upon these values exists, in the following we 
will assume A ~ 0.3. 

We use exact metric models, built by Occhionero et 
al. (1981,1983), for the description of spherical, compen- 
sated structures directly in the matter dominated era 
(MDE) . We investigate the metric distortion induced on 
the CMB radiation by the voids and check our approach 
by comparison with previous analyses of strongly non- 
linear highly approximated voids completely crossed by 
photons (Thompson & Vishniac 1987, hereafter TV). 
CMB complete crossing of voids was treated recently by 
some authors (Panek 1992, Arnau et al. 1993). We find 
that the metric CMB distortion in the case of complete 
crossing, Rees-Sciama effect (1968, hereafter RS), is neg- 
ligible with respect to partial crossing, Sachs- Wolfe ef- 
fect (1967, hereafter SW), for voids intersecting the last 
scattering surface (LSS), the width of which is modelled 
by making use of exact decoupling physics. Therefore, 
for treating accurately this case, we evaluate the CMB 
adiabatic distortion before decoupling due to the void. 
This is achieved by adapting the formalism developed 
by Hu and Sugiyama (1995, hereafter HS) to our case: 
it is shown that on scales in which SW effect is already 
comparable with the detected CMB anisotropics, the 
adiabatic effect is significatively reduced by Silk damp- 
ing and then may be neglected. These results have been 
obtained by recovering, both for the SW and the adi- 
abatic case, the angular dependence of the CMB tem- 
perature, which is therefore a function of the void size, 
of its exact density profile and of the particular position 
of the void with respect to the LSS; our model responds 
to all these items. 

We anticipate that the largest void on the LSS that 



does not violate the COBE observations has a comoving 
size of ~ 25 /i _1 Mpc, i.e. physical size 100 h _1 Mpc by 
today due to over comoving growth. 

The content of this paper is as follows. In Sec. II 
we construct null geodesies . In Sec. Ill we analyse the 
CMB distortion induced by complete crossing of a void, 
comparing with previous results based on highly approx- 
imated models. Finally, in Sec. IV we consider the obser- 
vational consequences of the voids lying on LSS. In the 
last section we draw the conclusions. 



2 GEODESICS 

The most natural theoretical enviroment to describe 
non-linear spherical large-scale perturbations was found 
by Tolman in 1934. We make the minimal assumption 
of pressureless fluid, according to the CDM paradigm. 
The fundamental length element is (c = 1) 



( dr Y 



T 2 (M) \dM ) 
- r 2 (M, t) (sin 2 6dcj> 2 + d6 2 ) , 



(5) 



where M, the gravitating mass, plays the role of a radial 
coordinate, and T(M), a generalization of the Lorentz 
7, is a conserved quantity in the pressureless case; it is 
given by 

dr\ 2 2GM 



I'M 
Also 

p(M,t) = 



1 



(6) 



(7) 



4nr 2 dr/dM ' 
The r(M) function characterizes the perturbation 
structure. The choice V = 1 reproduces the unperturbed 
Friedmannian case. For hyperbolic perturbations, we 
write 



T 2 (M) = i + r 2 h (M), 

and the solution to Eq. (^) takes the form 
GM 



r{M,t) 



t = 



Y\{M) 
GM 



(coshr/ — 1) 
(sinhr; — 77) 



(8) 

(9) 
(10) 



n(M) 

The model that we use here for Y\(M), Occhionero et 
al. (1981,1983), is suitable for description of a wide class 
of compensated matter structures, and has the following 
analytical form: 
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Tl{ X ) = B(±) X ^ n { X ) , *=ii , (11) 
/ 3M. ^ V3 



L s = 2GM„ , L 



\47rpcx3,o 



(12) 



7n (z) = e X p[-(r E (l/n)/n)"x W3 ] (13) 
(r^ = Euler Function). 

The model contains a tunable scale M, that will be typ- 
ically 10 15 Mo for a large scale void. The BL S /L factor is 
a measure of the perturbation strength. The exponential 
in Eq.([l3]) sets to zero the perturbation for M 3> M»; 
7n(a;) becomes a sharper and sharper stepwise function 
by increasing n, while simultaneously the compensat- 
ing shell becomes narrower and denser; the coefficient 
in the exponent is dictated by the request of measur- 
ing the energy associated with the perturbation. With 
this model for V+(M), the dynamics is obtained solving 
Eqs. (|),©. 

We have now a model of spacetime; the next step is 
to integrate the null geodesies. We place the observer on 
the direction <f> = 0, 9 — tt/2 from the void's center. A 
point on a null geodesic is marked by an affine parameter 
I, and the wave vector is 

(K t ,K M ,K e =0,K*)=±(t,M,0 = l4>) . (14) 
The geodesic equations are 



dK l ^ 1 dr d r , M n2 _ dr /„^\2 

di r 2 dM dtdM v > r dt^ > 



(15) 



dK 
dl 



KdMJ { > r dM v ' 



2 (dr_\ 1 &2r K t K M 

KdMJ dtdM 
( dry 1 d 2 r ,„ M - 
KdMJ dM 2 



dK* f dr^ m 9r , 

dl r \dM dt 

Of these, Eq.(^) integrates directly into 



(16) 
(17) 

(18) 



where the constant k is fixed by the initial conditions. 
We now integrate numerically Eqs.(|l^)([l6]) and ( p^ ) us- 
ing the normalzation 



K a K a = 



(19) 



as a check for the numerical accuracy. We set the "ini- 
tial" conditions at the observation point Pq and we inte- 



grate backwards. After choosing (K t )o = —1, Eq. Q1S|) 
at Pq reads 



(K M ) 



-r(—Y* 

KdMJ 
{3M\ 



* , (20) 
V r / 

where in the last equality, obtained by imposing that 
the observation point lies in the Friedmann region, a is 
the angle that the geodesic under investigation forms at 
Po with the direction of the void's center: 



k — rosina 



(21) 



With the relations K 4, = (l/r) sina and K B = 0, the 
initial conditions are complete. A similar approach to 
the construction of null geodesies on spherical structures 
was used firstly by Panek in 1992. 

Finally, we know (see e. g. Anile & Motta 1976) 
that the CMB temperature measured by the observer 
is related to its value seen at any point x along the 
geodesic specified by a through the relation 

To ' 



(22) 



where u M is the four-velocity of detection instrument. In 
the most interesting situation of a void intersecting the 
LSS, a peculiar velocity field causing a Doppler contri- 
bution exists at emission due to acoustic oscillations of 
the photon-baryon plasma: this contribution, evaluated 
in Sec. (IV) together with the dominant or comparable 
adiabatic effect, is small with respect to SW one because 
of Silk damping. 

3 VOIDS IN FRONT OF THE LSS 

As we already observed, a void, whether compensated or 
not, shows an over-comoving growth in an ordinary flat 
Friedmannian background. Then, in a first approxima- 
tion, a void can be schematically represented as a region 
of space in which the scale factor a(t) grows faster than 
t 2//3 ; more precisely, Bertschinger (1983,1985) found t 4 / 5 
for compensated voids in the highly non-linear, asymp- 
totic regime; our structures, being mildly non-linear, ap- 
proximatively obey this law. Such a dynamics generates 
a cooling (redshift) in the CMB temperature (in addi- 
tion, of course, to the ordinary one due to the unper- 
turbed cosmic expansion). 
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Figure 1. RS effect: void between us and the LSS; complete crossing of the void by the CMB photons. The non-linear density 
profile of panel a induces a very characteristic angular dependence on CMB temperature (panel b). In panels c and d dependence 
on R (~ R [i ) of the central redshift and of the peripheral blueshift. 
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Figure 2. RS effect for different values of n (different sharpness of the density profiles) and for interesting values of the present 
radii. The values reported for < ST/T >, practically indistinguishable for n=l and 2, are obtained by simulating COBE's 
angular resolution and appear much below the observed anisotropies. 
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On the other hand, the radiation is heated when 
crossing the void walls, as it is discussed below. 

The case to be discussed here concerns a void at 
z v — 900, completely crossed by a CMB photon. This 
case was analyzed in TV (firstly) and in recent works 
(Panek 1992, Arnau et al. 1993, Martinez-Gonzales et 
al. 1990). We leave to the next section the more impor- 
tant case of voids intersecting the LSS. TV found an 
analytical form for the angular dependence of the CMB 
temperature by adopting a model for the metric inside 
a compensated void which, albeit very schematic, con- 
tains all the essential aspects of the problem. As it is 
shown below, our result remarkably agrees with theirs. 
Their void is simply a spherical Friedmannian region 
of space in which the scale factor grows in time faster 
than in a fl = 1 universe, precisely as i 4//5 . With an in- 
genious double Lorentz transformation, they find that 
photons entering and leaving the void are boosted by 
Doppler effect, simply because they are seen by frames 
moving towards them; on the contrary, during the void's 
crossing, they cool with respect to the flat Friedmannian 
case, because the space growth is faster. As the second 
effect is strongly dependent on the time spent by pho- 
tons inside the void, while the first is almost fixed, the 
blueshift becomes dominant for strongly off axis trajec- 
tories. 

Thompson & Vishniac's analytical result is: 

ST _ f R/(l + z) ,x ' 
T 



H- 



- - cos 2 Vc) , (23) 



where ip a is the angle between the perpendicular to the 
shell and photon's direction at the leaving point; its ex- 
pression as a function of a is 



+COS 



do sin a 

r/(i + zy 



+ cos a\ 1 — 



(24) 



where d v is the comoving distance of the void's center. 

In Fig.(jl|) we report our numerical analysis of a 
typical void along the line-of-sight (R = 30 /i _1 Mpc 
, z v = 900). Panel (a) shows the density profile as a 
function of the radial coordinate M normalized to M„ 
the latter being calculated through the value (~ -R/4) of 
the comoving radius at the time of interaction with the 
CMB. The central underdensity, obtained for B = 2200 
in Eq.rtlll), and the large compensating shell are well 



visible. Panel (b) reports the CMB angular temperature 
change for complete crossing; the curve presents all the 
discussed features and agrees with the result of Eq. (^3|) . 
A completely crossed void appears in the CMB sky as a 
cold circular spot sorrounded by a hot ring, because in 
off axis trajectories the hole contribution is suppressed. 
The value predicted by Eq. (|2^) for central trajectories 
is — 10 -4 ; our result is slightly less because our voids are 
not completely empty. The maximal redshift to maximal 
blueshift ratio is —2.4 in Eq. (j^) as in our case. Finally, 
note the smallness of angular sizes, due to the fact that 
this void is at the edge of the observed universe; the 
temperature distortion induced by such a structure is 
much below the experimental limit (~ 10~ 5 ) because 
the present instrumentation has a finite angular width 
(3° Gaussian beam for COBE). Panel (c) and (d) show 
the dependence of maximal redshift and blueshift on 
the void's radius. The solid line is numerical, while the 
dashed one, pratically indistinguishable from the first, 
is {R/{l + z\ 2 /H- 1 f. 

Now, let us come to the evaluation of the CMB de- 
pendence on the shell sharpness. As we stated above, 
the latter is regulated by the parameter n. We have 
confronted the CMB temperature changes induced by 
different n's, and for voids of different sizes in the inter- 
esting range. The results are shown in Fig.(^) and reveal 
a very weak dependence on n. The same happens dur- 
ing void crossing. In this case (the most interesting) , for 
any value of the parameter BL S /L in Eq. (|ll|), large n 
makes 7 n in Eq.(|l3|) approach the step function; in turn 
that yields the limiting sharpness for the shell's density 
profile. As we found that the variation of n does not 
change significatively the CMB perturbation (Eq. (^) 
below), we will neglect in the following this effect, and 
work with the simplest case n = 1. Also, Fig. (Q) shows 
< 5T/T >, defined as a COBE-like Gaussian-averaged 
signal: 

<^>=7V J ^(a)exp[-Q 2 / 2 7 2 ]27rsinadQ , (25) 

where 7 = 3°, COBE's beam size, and N is the nor- 
malization constant; in the following, we will refer to 
the COBE Gaussian as W(a). It is evident that, as ex- 
pected, the values reported in Fig.(^), of the order of 
10 -9 for the largest R, are extremely small compared 
with the observed anisotropies. 
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Figure 3. Temperature changes followed as functions of z along a central trajectory in voids with present physical radii R =30, 
50, 80, 100 ft Mpc sitting exactly on the LSS. In panel a the photons come from infinity with vanishing ST/T; realistically, in 
panel b they leave either at z = 1100 or at z = 900; it is then seen that their frequency shift is cither toward the red or toward 
the blue depending on the departing position. It must be remarked that the order of magnitude of this metric induced change 
is nearly 10 3 times its corrispondent value for complete crossing. Finally, panel c reports the underlying density profile. 
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We underline that this extreme smallness derives 
from the complete crossing of the void: in fact, the 
only observationally interesting voids are those inter- 
secting the LSS, to be analyzed below. Furthermore, it 
must also be underlined that voids with z <C 1000 will 
produce an even smaller temperature perturbation, be- 
cause cosmic expansion, not compensated by the over- 
comoving growth, dilutes them. 

4 VOIDS ON THE LSS 

Fig- (Q) shows the main new result of our work. In the 
upper and in the central panels we report the temporal 
evolution of the photon temperature during its crossing 
of the void: the latter is shown in the lower panel; differ- 
ent void radii of interest for the observations are shown. 
The abscissa is z — zlss, the latter being the redshift of 
LSS (~ 1067); in the four cases shown the void's center 
lies exactly on the LSS. In the upper panel the photons 
come from infinity with ST/T = and experience firstly 
the blueshift from crossing the first shell, secondly the 
strong central redshift generated by the over comov- 
ing expansion, and finally the blueshift from crossing 
the second shell: globally this yields the small result- 
ing ST/T (RS effect) hidden on the left hand border of 
curves in Fig.(|^). This is a very characteristic image of 
a void crossing, because we have found it being nearly 
the same in form for all the trajectories and astrophys- 
ically interesting void sizes. The central panel applies 
to realistic paths starting from either edge of the LSS 
(z ~ 900 and 1100) and displays the total blueshifts 
and redshifts that can be read on the left hand side. 
The order of magnitude of the temperature change dur- 
ing crossing is predicted by considering the propagation 
of a photon in a region of space expanding as t 4 ' 5 during 
a time St ~ R; the result is 



ST 
T 



R 



R 



5(1 + z)- 2 H- 1 - 20 H- 1 ' ^ 
where overcomoving growth is taken into account and 
is calculated at interaction epoch z ~ 1000. As 
may be easily seen from the figure, the predicted value is 
strictly respected; this argument was inferred firstly by 
Vadas (1995). For the reason that ST/T is proportional 
to the time spent in the void, for non-central trajecto- 
ries 7? must be replaced with Ry/l — (a/an) 2 , just the 
void's length on a direction characterized by the angle 



a defined in Sec. (II), the maximum value of which is 
olr. An accurate analytic approximation of the results 
for all trajectories is 



ST _ R y/l-[a/a R y 



T 

where 



20 H- 



1.25 x exp 



(— V 



-(27) 



(28) 



y/1 - (a/a R ) 2 R/4 
and r is the photon's coordinate. We see that the sharp- 
ness is already very strong in the present case n = 1; 
increasing n up to infinity simply brings the curves in 
Fig. (til) to their limiting case 



ST _ 
T ~ 
where 



R y/1 -(a/an) 2 



1.25x0(|a;| - 1) 



(29) 



20H- 1 

is the Heaviside function. We also see that an 
axial trajectory presents a typical ST/T three orders of 
magnitude above the central one of Fig.([j]) (panel (b)), 
obtained for complete crossing and hidden on the left 
hand border of curves of the RS effect, displayed in the 
upper panel of Fig.(^). The reason why this strong con- 
tribution does not appear in Eq. ( p3| ) is that for com- 
plete crossing the effects of the shells and the central 
hole balance out, as it is evident from Fig.Q, leaving 
place to terms with higher power of R/H -1 . Therefore, 
if a void intersects the LSS, the ensuing ST/T may be 
very different depending on the particular LSS position 
with respect to the center of the void. In any case, the 
amplitude will be much larger than in the case of com- 
plete crossing. Let us make a simple calculation: an use- 
ful approximation to the total temperature distortion 
induced by a void of angular size an averaged over a 
beam of size 7 is 



and taking ST /T\ a= o ~ 10 -2 as a typical value from 
Fig.(g) we obtain for R = 100 /i _1 Mpc (R = 25 /i _1 Mpc 
at decoupling) and 7 = 3° 

<^>^6-10- 5 . (31) 

This value, barely above the present observations, indi- 
cates that such scales are very close to the maximum 
permitted to primordial voids from CMB isotropy. This 
motivates the accurate investigation below. 
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Figure 4. Purely adiabatic and Dopplcr ST/T due to a compensated void in the cosmic fluid, both in the approximate (solid 
curve) and exact (dashed curve) case. On the relevant scales, this perturbation is well below the SW effect. 
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To treat accurately the case of voids intersecting 
the LSS, we must take care of all the physical processes 
at work. In fact the SW effect on the photons decoupled 
inside a void is not the only perturbation to the CMB; 
the corresponding inhomogeneity of the photon-baryon 
plasma naturally induces adiabatic and Doppler tem- 
perature fluctuations before decoupling. We anticipate 
that such contributions are reduced by Silk damping 
on scales < 25 /i -1 Mpc in which the SW effect is al- 
ready comparable with the observed CMB anisotropy. 
To see this, we have used the formalism developed by 
HS; it allows to obtain the CMB fluctuations at decou- 
pling caused by a perturbation characterized by a New- 
tonian potential <f) (our perturbations are not properly 
in the linear regime, but this approximation allow us 
to estimate correctly the effect). In the HS treatment, 
the Boltzmann, Euler and continuity equations for the 
photon-baryon plasma, accounting for Thomson diffu- 
sion together with gravitational and pressure terms, are 
analytically solved under the tight-coupling approxima- 
tion. In our treatment, in order to estimate only the 
adiabatic and Doppler effects, we have excluded the SW 
terms from HS formulas; we performed the calculus ex- 
actly and with the following simplifications: a) constant 
gravitational potentials, which is nearly true in the lin- 
ear regime and approximatively in our mildly non-linear 
case; b) the quantity 3pt,/4p 7 , that goes from to .1 from 
the remote past to decoupling in models with Qb — .01, 
is neglected (equivalent to have a constant sound speed) ; 
c) istantaneous decoupling, and the Silk Damping scale 
taken at this time. With these hypotheses in the HS 
treatment, the amplitude of Fourier k-mode of temper- 
ature fluctuation at decoupling (_d) (we report only the 
monopole term for simplicity) is related to the potential 
by the simple and intuitive relation 



'ST-- 

v T / k.D 



, r, ! 1 - icosfcc^o) e - (fc/fcB)2 , (32) 



where \ D = 2n/k D = 3(ni,/i 2 )" 1/2 )(n o /i 2 )- 1/4 Mpc is 
the Silk damping scale at decoupling (again from HS 
formulas), and r\ is the conformal time. On scales larger 
than the sound horizon c s t]d, Eq.(|32"|) reproduces the 
known behaviour ST/T oc <j> up to a factor of the or- 
der of unity; on scales smaller than 1/ko Silk damping 
forces to zero the inhomogeneities of the photon-baryon 
plasma. At intermediate scales, corresponding to .2°-^l° 



in the sky, an oscillatory behaviour occurs due to the 
acoustic oscillations; it will be retained in the form of 
Doppler peaks in the correlation function of 8T/T. As in 
HS, one can antitransform (ST/T)k to obtain the contri- 
bution to ST/T from the adiabatic effect as a function 
of the angle a. For the potential source we have used 
the simple model 



dp, , 
— (r) 

P 



= ~ 5 > R^ 1 



Sp 



(<•) 



SP, s 

-(r) 



(1 + a) 3 - 
0, ^>l+a, 



1 < ^ < 



(33) 



i.e. a central cavity of depth 1 — 5 surrounded by a com- 
pensating thin shell AR — a ■ R. The result, together 
with the exact calculus, is shown in Fig.(^) for the in- 
teresting values (1 + zd)' 2 R = 25 /i _1 Mpc, a=.5, 8 = .7 
and the typical set fit = .06, Qo = lj h = .5 for a CDM 
universe; the signal that COBE would detect for this 
perturbation is also reported and proves to be below 
the SW contribution, as we have anticipated. 

Therefore, the picture of a primordial void sitting 
on the LSS which emerges very clearly from this analysis 
is that of a strong spot of metric origin and of either 
color on very small angular scales. 

For the above arguments, we are also allowed to use 
the standard homogeneous model for the width of LSS. 
By taking into account the dependence of the fractional 
ionization on redshift z, and its consequences upon op- 
tical depth for photons, we are led to the probability 
function for a photon to be last scattered at redshift z 
(Jones and Wyse, 1985) 



P(z) = 



5.26 



exp 



-0.37 



1000 



(34) 



1000 U000, 

that is well approximated by a Gaussian centered in 
Zmax ^ 1067 and with a ~ 90. 

Now, let's use current observational data on the 
CMB isotropy to estabilish constraints on voids at de- 
coupling. As we stated above, the qualitative form of 
the curves in Fig. (j^) is the same regardless of the 
voids' radii. The COBE observed anisotropies distribu- 
tion have an amplitude a ~ 3 ■ 10 -5 at small angular 
scales (Bennett et al. 1994): therefore we will assume 
that a void is excluded by current observational data if 
its < ST/T > is 3cr ~ 10" 4 . 
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Figure 5. Observational constraints on the parameters p,Rm- the long-dashed line refers to the request that a reasonable 
fraction of space (50% ) is filled by primordial voids; the short-dashed line is based on small hot spots statistics; the solid line 
is the present result: the number of voids inducing a visible CMB perturbation must be smaller than one. The shaded region 
gives the acceptable parameter space for primordial void distribution in the whole observed universe. 
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We have integrated the CMB fluctuation for many 
void radii in the relevant range, to get < ST/T > as a 
function of R: 

ST 



< f (R) >- 



2tt sin adadz^(R, a, z) W(a)P{z) .(35) 



We can immediately obtain an upper limit for R as fol- 
lows. At radial distance (1 + z)' 2 r pea k = R inside the 
void, the CMB distortion, whether red or blue, reaches 
its maximal value before going to zero (see Pig. (^J)) . If 
T P eak coincides with the peak in P(z) (at zlss ~ 1067), 
the CMB distortion detected by an instrument with 
Gaussian window function with width 7 (expressed in 
degrees) is found to be 



< 



5T 
T 



10" 



m ■ 



(36) 



where the numerical factor contains a correction 1/3 
due to the LSS width for R near to H~ l . An R comes 
from the SW effect, while the remaining R 2 is due to 
the beam's angular width. Then, COBE-like experi- 
ments (7 = 3°) say that no voids with R — .97 H~ x ~ 
100 /i _1 Mpc (25 at decoupling) or more exist on the 
LSS. We point out that the same result without consid- 
ering the LSS width would be R — .63 

Now that we have the largest void passing the CMB 
bounds, we pass to evaluate the statistical consequences 
on distribution ([!]). The position of the LSS peak with 
respect to the void center corresponds to a value of 
x defined in Eq.(^8|); for any x, there will be a ra- 
dius Rmax{x) for which the CMB distortion detected 
by COBE would be nearly 3cr; this is obtained numer- 
ically from Eq.(p5|), putting the LSS peak at r = x R 
from the void's center. The number of voids with co- 
moving radius R is \dNs /dR\dR, and the fraction of 
them having the LSS peak at distance between r and 
r + dr is 



dN(R,x) = 



dN E 



dR 



dR- 



3dr 
2H" 1 



dN B 



dR 



RdR 



3 dx 
2H- 1 



.(37) 



The total number of such voids with R > R ma x must 
be less than 1; this is equivalent to the condition 



dx 



dR 



„0O 



dN B 



dR 



3R 
2H- 1 



< 1 



(38) 



By using Eq.(gJ) we obtain the following upper limit re- 
lation between Rm and p: 

/V- ( ~ ^ I Rma.ix) 1 -*^^) . (39) 



2Ho 



The result is shown in Fig.(g) as a solid line; values 
above the long-dashed line ensures that at least half 
of the space is presently filled by large scale voids 
(R > 10 /i -1 Mpc ); values below the short-dashed line 
ensure that the random noise produced by an overposi- 
tion of the many small voids on the CMB does not break 
the CMB isotropy ( Amendola & Occhionero 1993) . The 
dashed area contains roughly the set of (p, Rm) values 

6 < P < 13 , 30 /i _1 Mpc < R M < 130 /i _1 Mpc , (40) 

that is the parameter space for which the voids respect 
all the CMB constraints worked out so far, and pro- 
vide large scale power matching current observations of 
galaxy clustering. 



5 CONCLUSIONS 

We have built a formalism that allows us to check 
the plausibility of the cosmogonic role of primordial 
voids; here, strictly speaking, by "primordial" we mean 
"present already at decoupling": however the obvious 
suggestion is that the voids originate much earlier, pre- 
cisely in an inflationary first order phase transition. The 
constraints from the observations lie entirely in the com- 
patibility of the voids' theoretical imprints on the CMB 
with the observed anisotropies at the relevant scales and 
translate into defining an allowed region of the param- 
eter space for the void distribution. 

In order to compute photon trajectories and tem- 
perature perturbations, we describe a void by an appro- 
priate analytic general relativistic model of a spherical 
perturbation embedded in a flat Friedmann background. 

Firstly we remind from the literature the case of 
a void lying between us and the last scattering surface 
(LSS): this implies the complete photon crossing known 
as Rees-Sciama effect. It is known that this induces a 
<5T/T~ -(4/15)( J R/(l + 2)- 2 ^ -1 ) 3 : we have found that 
this results from a partial cancelation of the redshift 
occurring during the void crossing, namely the Sachs- 
Wolfe effect, with the blueshifts occurring at either shell 
crossing. Each of these effects, \5T/T\ ~ R/20H' 1 , is 
three orders of magnitude larger than the total effect. 

In the present paper, we have therefore considered 
the case of voids lying on the LSS and we have eval- 
uated the global effect for photons decoupled inside 
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the void, comprising the Sachs- Wolfe contribution men- 
tioned above, the adiabatic and Doppler effects due to 
the acoustical oscillations of the photon-baryon plasma; 
the latter has been found small with respect to the 
former because Silk damping is active on the small 
scales involved. For completeness, we have used an ex- 
act model for the LSS profile. Therefore, we are led to 
a very clear picture for the CMB impact of primordial 
voids: they are strong spots of metrical origin on very 
small angular scale. 

A first observational constraint is obtained by re- 
quiring that the number of astrophysically interesting 
voids intersecting the LSS and capable of generating 
a distortion above the present observational limits on 
small angular scales, is smaller than unity. In practice 
this yields the available region of the parameter space 
for the void distribution, which we choose in the conve- 
nient form of a power law. 

Where do we go from here? We expect a lot of power 
on the very small angular scales, £ ~ 1000 in the canoni- 
cal expansion in Legendre polynomials of the CMB cor- 
relation function. Power on larger scales, £ < 100, comes 
from the ordinary zero point fluctuations of the inflaton 
as well as from the variation of the number of voids in 
different regions of the sky, while in our case an un- 
expected strong correlation of metric origin arises on 
scales where Silk damping erases everything else. We 
plan therefore to build realistic sky maps by taking into 
account all the voids on the LSS in order to recognize 
whether the imprint of primordial voids will be actually 
present in the data of future high resolution experiments 
like COBRAS - SAMBA. 

We thank David Wands and Michael Sahzin for use- 
ful comments. 
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